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We revisit the classical instability of closed strings in the Abelian Higgs model. The 

instability is expressed by a vanishing energy as the torus-like configuration shrinks to zero 

size. We treat the radius of the torus as a collective coordinate and demonstrate that a 

quantum mechanical treatment of this coordinate leads to a stabilization of the closed string 

C*~) at small radii. 
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I. INTRODUCTION AND MOTIVATION 

The study of string or vortex type configurations in the Standard Model and related theories 
has attracted much attention in the recent past [THS]. Because these configurations may have 
cosmological implications and to distinguish them from the fundamental objects in string theory, 
they are commonly called cosmic strings. The most imperative question regarding the dynamics 
of a cosmic string is its stability. Studies of this question have mainly focused on infinitely long 
and axially symmetric string configurations [&M22j . Unless these configurations carry (topological) 
charges, as is the case, for example, for the Abelian Nielsen-Olesen string [6], they are classically 
unstable but can eventually be stabilized by quantum effects. One such possibility would be to 
bind (a large number of) fermions at the center of the string |21|l22j. Other stabilizing mechanisms 

qq involve the coupling of external fields or currents leading to superconducting strings [23]. 

r ■ Strings at the boundary of regions with different vacuum expectation values are expected to 

build a network, so that they have the potential to form cosmic string loops |23} I24j . As a con- 
sequence, it is conceivable that closed strings or a network thereof existed at early times of the 
universe when the vacuum manifold was not simply connected. Cosmologically motivated studies 
on networks of closed strings suggested their stability (against graviational radiation) for radii of 
the order of a a small fraction ~ ^(0.01 . . . 0.1) of the horizon distance. These studies started from 
analytic investigations |25H27j but were followed by intensive numerical simulations, cf. the recent 

^ article [28] for a comprehensive discussion of these activities over the past two decades. 

Eventually such closed strings would relax to circular (torus-like) configurations, so-called vor- 
tons [21] . Typically such vortons are stabilized by coupling to (external) currents [23] [29H33] . The 
existence of vortons in a super-conducting scenario was supported by numerical studies some time 
ago [Ml. 

It is also interesting, however, to see whether such closed strings are self-stabilizing in the sense 
that the dynamics of the underlying Lagrangian prevents them from contracting and eventually 
disappearing. We are particularly interested in quantum self-interactions of a single closed string. 
We expect that these quantum effects become significant for strings with radii similar to the 
Compton wave-lengths of the fields within the Lagrangian. In such a scenario the string would lose 
energy by emitting fluctuations of its own field rather than by graviational radiation. 

Closed strings are not subject to any non-trivial boundary condition at spatial infinity and are 
therefore not stabilized topologically. On the contrary, simple dimensional arguments show that 
closed strings can reduce their energy by shrinking. Hence they are not stabilized dynamically 
either, at least on the classical level. The naive conclusion would be that even if a closed string 
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had existed at some earlier time it would have decayed by emitting radiation and simultaneously 
shrinking to zero size. However, this argumentation cannot be correct at the quantum level because 
narrow configurations localize the particle in position space which, by the uncertainty principle, 
implies that large momenta emerge and prevent the total energy from vanishing. 

We will consider this question for the flux tube in scalar electrodynamics and study the analog 
of the Nielsen- Olesen string on a torus of radius R. In this case a classical instability occurs as 
R — > 0. The study of a toroidal string configuration is significantly more complicated than that 
of an infinitely long and axially symmetric string because the field equations cannot be simplified 
to ordinary differential equations of a single variable. This hampers the application of spectral 
methods, which have previously been used successfully to compute quantum effects for localized 
field configurations [35]. We will therefore adopt a different technique that has proven success- 
ful in describing similar quantum mechanical scenarios. It is based on the quantum mechanical 
uncertainty relation 

ApAx > - (1) 

and assumes that for the ground state the momentum p is saturated by Ap and that Ax is a good 
approximation for the spatial extension, x. Minimizing the Hamiltonian that is constructed from 
the replacement p — > h /{2x) yields the correct ground state energy and extension for both the 
hydrogen atom (which is also classically unstable) and the harmonic oscillator. 1 

In order to apply this approach, the appropriate momentum (operator) must be identified. We 
will do this by introducing a time-dependent collective coordinate for the extension of the field 
configuration in which the classical instability occurs. The conjugate momentum is then extracted 
from the Lagrangian and the system is quantized by imposing the canonical commutation relation 
between the collective coordinate and its conjugate momentum. This so-called collective coordinate 
quantization reasonably reproduces quantum properties of localized field configurations in various 
examples [36]. Rather than exploring the full quantum theory, this formulation concentrates on 
the physically important modes. In the case of a collapsing cosmic string, the most relevant mode 
is obviously the radius R of the torus on which the string lives. 

The paper is organized as follows. First, we complete this introduction by describing the model 
of scalar electrodynamics that we use. In section II we briefly review the Nielsen-Olesen string 
in the axially symmetric framework. In section III we provide a detailed description of toroidal 
coordinates and explain the string configuration that respects this symmetry. In particular, we 
derive consistent boundary conditions. As expected from dimensional considerations, we find that 
the classical toroidal string is unstable against shrinking to zero size. We therefore introduce 
a collective coordinate for this mode, which we quantize canonically (section IV). We estimate 
the corresponding contribution to the energy from Heisenberg's uncertainty principle and derive 
a quantum energy functional whose minimum we seek. Since the system cannot be decomposed 
into one-dimensional problems, the numerical investigations are challenging. We relegate their 
description to appendices. Section V summarizes and concludes our studies. 

We consider scalar electrodynamics with the Lagrangian 

if = -± Fflu F^ + \D fl <Z>\ 2 -V($), (2) 



1 For the hydrogen atom one needs to replace Ax — > 2r, where the radial coordinate r is the distance from the center 
of the Coulomb potential. 



3 



where the covariant derivative and scalar potential are 

D li = d li + ieA li and V($) = m 2 |3>| 2 + A |$| 4 , (3) 

respectively. To induce spontaneous symmetry breaking, we take m 2 < and observe the classical 
vacuum expectation value (vev) 



i*i> = y^». w 

Notice that $ and v have dimension of mass in (3 + 1) dimensions, while A and e are dimensionless. 
The (tree-level) masses of the Higgs and gauge fields can be identified by considering fluctuations 
about the vev in the Lagrangian in eq. Q, 

m| = 4Xv 2 , m 2 A = 2e 2 v 2 . (5) 

For later reference, we briefly discuss the equation of motion for the temporal component of the 
gauge field (Gaufi' law), 

( - A + 2e 2 |$| 2 ) Aq = ie ($*d $ - <S>d <Z>*) + V • A , (6) 

where the dot denotes a time derivative. Clearly, this equation contains no time derivative for Aq, 
and is therefore a constraint on Aq rather than a true dynamical equation. Furthermore, it is an 
elliptic differential equation and as such has the unique solution Ao = (given proper boundary 
conditions at spatial infinity) if the right hand side of eq. ^ vanishes. This happens if (i) the 
phase of the Higgs field is time-independent and (ii) the gauge field A is either time-independent 
(A = 0) or transverse (V • A = 0). It is, of course, possible to enforce the transversality of 
A (or the Weyl condition Aq = 0) by a choice of gauge, even for time-dependent fields. The 
respective gauge transformations will, however, be time-dependent and highly non-local, leading 
to very cumbersome contributions to the phase of the Higgs field and the kinetic energy of both 
A and 3>. Since the kinetic energy is essential for the potential quantum stabilization of a closed 



string (cf . section IV ) , it is very inconvenient to enforce the Coulomb or Weyl condition in order to 
resolve Gau8' law. Instead, we have to implement eq. Q and allow for Aq ^ when introducing 
the time-dependent field configurations necessary for a quantum mechanical treatment. 



II. INFINITELY LONG NIELSEN OLESEN STRING 

Before discussing our toroidal string configuration in more detail, let us briefly review the well- 
known Nielsen-Olesen (NO) string [6j which corresponds to the limit that the (outer) torus radius 
becomes infinitely large. This string is a static, infinitely long vortex configuration stretched along 
the z-axis, in which the Higgs field <3? has winding number n, and the magnetic field points along 
the string axis, B ~ e z : 

A = g(r)-e v , ^ = vf(r)e+^. (7) 

er 

(We are using polar coordinates (r, ip) in the xy-plane) . Notice that V x A = and Aq = for 
this ansatz, which is consistent with Gaufi' law because the Higgs field $ is time-independent. The 
requirements of finite energy and continuity at the origin then lead to boundary conditions on the 
profile functions / and g, 

r -> : f(r) -)■ 0, g(r) ^0 

(8) 

r — > oo : f(r) — > 1 , g(r) — > 1 
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The NO string owes its stability to the fact that the Higgs field at spatial infinity defines a map 
S 1 i — y U(l), which in view of 7Ti(£7(l)) = 7L cannot be continuously deformed into the vacuum. 

We note that our model, as given by eq. Q, initially has three parameters, the Higgs vev v and 
the two coupling constants A and e (or, equivalently, the two masses in eq. ([5])). After rescaling all 
dimensionful quantities in units of some appropriate scale, we expect all observables to be functions 
of the two couplings A and e independently. However, this is not always true: In the NO string, 
for instance, the profile functions only depend on the ratio 



/3=4 = ^ <») 



once the fundamental scale is chosen as the mass of the gauge field, (and the winding number 
is fixed to n = 1). This is because the gauge coupling e scales out in an overall prefactor for the 
(dimensionless) energy per unit length, 



E tt f , ( n 



m l A L e 2 J \ p 2 



o 



,2 
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dpp\-,g'{pY+ f(p)* + -f(pf (l-g(p)Y +£(1 -/(,)')* , (10) 



P 2 
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where p = r is the scaled distance from the symmetry axis. The NO profiles f(p) and g(p) 
resulting from the minimization of this functional will only depend on the ratio /3, because any 
change of couplings that leaves j3 invariant can at most produce an overall prefactor to the functional 
(10), which is irrelevant for the minimization. The value of the (minimal) energy per unit length 
will then scale trivially with e at fixed /3. In figure [TJ we present the numerical solutions for the 
NO profiles in the case j3 = 1, and the value of the minimal energy per unit length as a function 
of p at fixed e 2 = 1/2. For our preferred values 

e 2 = - , A = - => wia = rn<& = v and /3 = 1 , (11) 

we find a minimal energy per unit length of 6.28 m 2 ^; this will provide an important numerical 
check for the torus configuration discussed in the next section. 



III. TOROIDAL COORDINATES AND STRINGS 

We construct a closed string from the NO configuration by first assuming that the latter has 
a finite (though very large) length L = 2irR and then identifying its two end surfaces. The core 
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of the string, along which the profile functions / and g vanish, thus becomes a circle of radius R. 
Moreover, the radial coordinate p, i.e. the independent variable of the NO string, then measures 
the distance from this circle, and surfaces of constant p ^ are tori. We will be particularly 
interested in the change of the profile functions as the core radius R decreases and becomes as 
small as the Compton wave-lengths of the particles in the model. Eventually R will be considered 
a variational parameter which will be determined by minimizing the (quantum) energy of the torus 
configuration. 

The NO string is characterized by a rotational symmetry about its core. This cannot be 
maintained once the string is closed, because the direction in which one moves away from the core 
does matter: If one moves away in the direction of spatial infinity, the fields have an infinite range 
to decay to their vacuum values, while there is only a finite distance when moving towards the 
center of the core circle. (And it is even not necessarily true that the fields assume their vacuum 
values at the center.) As a consequence, the lines of constant energy density must be denser on the 
inside of the core than on its outside, and since these situations are related by rotations about the 
core axis, there cannot be any axial symmetry. While the profile functions at spatial infinity are 
still subject to boundary conditions, they will result from the dynamics at the center of the circle. 
The loss of axial symmetry causes the profile functions to depend on more than one coordinate, 
which complicates matters significantly. The only remaining symmetry is that along the core circle. 



A. Geometry of the torus configuration 

The previous considerations strongly suggest the introduction of toroidal coordinates which we 
describe in this subsection. We choose to put the circle of the string's core in the xy-plane of a 
Cartesian coordinate system and center it at the origin. The corresponding toroidal coordinates 
read [57]: 

sinh r 

x = R cos ip 

cosh r — cos a 

sinhr 

y = R — sin ip 

cosh r — cos a 

^ sine . . 

2 = i?— . 12 

cosn r — cos a 

The coordinate ip £ [0, 2ir] is the usual polar angle describing the rotation of the xz-plane around 
the vertical z-axis. By construction, this rotation is a symmetry of our configuration and <p is a 
cyclic coordinate that describes the position along the string core. The coordinate r E [0, oo] is 
best understood from its iso-lines in the left panel of fig. [2] Smaller values of r correspond to 
larger circles 2 that spread out towards infinity on the outside, and approach the vertical axis on 
the inside, while the center of these circles gradually moves away from the torus core. The inverse 
r _1 is therefore (highly non-linearly) related to the distance from the string core (similarly to the 
radius p in the NO string), although it cannot be directly identified with this quantity. 3 Finally, 
the angle a G [0, 2tt] describes the different points on the r-coordinate lines. When viewed from 
the core, it plays the role of an azimuthal angle. It is therefore natural to give the Higgs field 
winding by taking a as its phase. Then the profiles are predominantly functions of r only, with 



2 These are indeed circles: (for (p = 0) [x — -Rcoth(r)] 2 + z 2 = i? 2 cosech 2 (r) is independent of a. 

3 As we will see below, the closest analog to the radius p in the NO string is the metric factor h(a, r) explained in 
eq. ( 15 1 below. 
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x/R x/R 

FIG. 2: (Color online) Coordinate lines of toroidal coordinates in the iry-plane (<p = and tp = %). Left 
panel: Lines of constant r (at tp = only, <£> = 7r is the mirror image about the z-axis). Right panel: Lines of 
constant a. Each circle is composed of two distinct coordinate lines, namely a £ [0, 7r/2] (above the foci) and 
(er + 7r) (below the foci). The remaining a range gives mirror images which are not displayed for simplicity. 



a cr-modulation mainly in the vicinity of the origin. Note that origin of the coordinate system 
corresponds to r = and a = tt while spatial infinity is approached as r — > and a — > 0. 

The above considerations lead to the following ansatz for the closed string field configuration [38] 

$ = vf(a,T)e in ' 7 , A = A(a,T)e a , (13) 

where e a is the tangential unit vector along the lines in the left panel of figure [2} To ensure finite 
energy, we must have |D3>| —> at large distances, which implies 

, . . n 1 n cosh r — cos a ... 

A((T,T)^-- r = . 14 

e h(a, t) e R 
Here, the gradient in toroidal coordinates gives rise to the metric factor 

h(a,r)^—^ . (15) 

cosh r — cos a 

In view of these asymptotics, we refine our ansatz to the form 

^) = vf{a ) r)e m \ A (x) = -f^e CT , (16) 



with suitable boundary conditions on the profiles discussed below. Since the configuration eq. ( 16 ) 
is time- independent, we can discard any temporal gauge field component and put Aq = to resolve 
Gaufi' law. Note, however, that the gauge field is no longer transverse, 



„ . n 1 d 
V • A 



e h(a, t) 3 da i 

As expected, the magnetic field 



g(a,r)h(a,r) (17) 



B = VxA = -^e, (18) 
e h(a,T) 



of the closed string points in em-direction, i.e. along the core of the string. 
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B. Classical energy and boundary conditions 

Since our closed string ansatz is static, its energy density is just minus its Lagrangian density. 
The classical energy becomes, upon integration over space, 



- 2 \ dr da sinhr " (d T g) 2 + Rr, \(d T f) 2 + n 2 f 2 (l - gf + (d a ff 
3 io J-n Rn L 



+ ^ 3 r ? 3 (l-/ 2 ) 2 |. (19) 

As in the case of the NO string we have scaled all dimensionful quantities by appropriate powers 
of the mass of the gauge boson, m^: 

R^m A R, ^ a>T ) = M^l = . (20) 

R cosh r — cos a 

Again, the gauge coupling e only appears as an overall prefactor, and the profiles only depend on 
the mass ratio (3. In addition we observe a non-trivial dependence on the radius R. 

To complete the discussion of the classical torus configuration, we need to specify the boundary 
conditions for the profile functions / and g. Since a is an angle variable and the profiles have 
to be single- valued, it is immediately clear that we must impose periodic boundary conditions in 
cr-direction, 

/(0,r) = /(2vr,r), g (0, r) = <?(2tt, t) . (21) 

At r — > oo, we are on the core circle where we expect the initial U(l) symmetry to be restored, 
i.e. the Higgs vev vanishes. To avoid singularities at r — > oo, the magnetic field and even the gauge 
potential itself must also vanish on the core line. This enforces Dirichlet conditions 

f(a,oo) = g(a, oo) = 0. (22) 

Finally, we have to find the conditions for the limit r — > 0. This corresponds to the largest 
circles in the left panel of fig. [2] and includes spatial infinity (a — > 0, 2ir) and the vertical z-axis 
(a ^ 0,27r). The first case obviously requires the vacuum configuration / = g = 1. Furthermore, 
when approaching the z-axis from any nearby point identical field configurations must be produced. 
Hence all field derivatives must vanish as r — > 0. This imposes Neumann boundary conditions for 
a ^ 0,2tt, 

d T f(a,0) = d T g(a,0) = 0, (23) 

Any solution of the field equations with finite total energy will automatically approach / = g = 1 
as <t — } 0, 2tt when r = 0. This restriction results from the terms with positive powers of r/ in 



eq. (19), which would diverge otherwise. 



C. Numerics and classical instability 

Application of Derrick's theorem to the covariant Lagrangian shows that no localized configuration 
is classically stable, since there is no coupling of negative mass dimension. We will reproduce this 
classical instability explicitly by numerically calculating the minimal energy for a given radius R. 
This will provide us with important information and techniques needed to subsequently investigate 
quantum-mechanical effects. 



s 



To minimize the energy in eq. (19) subject to the boundary conditions, eqs. (21)-(23), we 



discretize the coordinate space by a square lattice of (N s + 1) x (N t + 1) grid points, 

T = T +t , t = 0,...,N t 

a = s — , s = 0,...,N s . (24) 

For our numerical treatment r = and r = oo have been replaced by finite values < tq <C 1 and 
Too 3> 1, respectively. The discretization replaces the continuous profiles by indexed quantities 

f st = f(sA s ,T + tA t ), g st = g(sA s ,T + tA t ) (25) 

defined on the lattice sites. The grid spacings are At = (r^ — To)/N t and A s = 2ir/N s , respectively, 
which can be made arbitrarily small by increasing the lattice size. With the integrals replaced 



by Riemann sums, the minimization of the classical energy (19) turns into a high-dimensional 
discrete minimization problem for the variables f s t and g s t, which we solve by a combination of 
standard algorithms such as iterated overrelaxation and simulated annealing. Further details on 
the discretization and the minimization algorithms can be found in appendix [TJ 

Our model has a number of parameters, which can be partitioned into two classes: 

1. physical parameters: These are the gauge coupling e, the ratio of couplings (5 given in eq. Q, 
and the dimensionless core radius R. As explained earlier, the gauge coupling only appears 
as an overall prefactor in the energy once dimensionless variables have been introduced. We 
will therefore adopt e 2 = 1/2 as a standard choice and study the energy as a function of (3 
and R. 

2. numerical parameters: These are related to the discretization and should become immaterial 
in the continuum limit N s — > oo and Nf — > oo. In addition to the grid size the numerical 
treatment also contains the cutoffs To and for the r-parameter. In practice, the cutoffs 
are replaced by 

2 2 

e T °° + 1 e T ° — 1 

where Rd m i n / mSLX are the distances of the r-isoline from the torus core, measured on the 
x-axis inside (cZ m in) or outside (d max ), cf. fig. [2j Hence the continuum limit also requires 
dmin -> and d max ->■ oo. 

First, we verify that the numerical (discretization) parameters become irrelevant in the contin- 
uum limit. This exercise also suggests appropriate data for the numerical parameters to be used 
in later computations. To this end, we hold all of them fixed at reasonable values 

N s = N t = 300, d min = HT 3 , (f max = 20.0 (27) 

and vary one at the time. The physical parameters are fixed to R = 10 and /3 = 1 in all cases, so 
that distances are all measured in units of the Compton wavelength of both the Higgs or the gauge 
boson (which are equal). 

Tab. [I] shows the results. The largest dependence is on the short distance cutoff d m \ n and the 
lattice size N s = Nt, but only if these are taken too far from the continuum limit. For conservative 
values d m ; n < 0.01 and N s = N t > 150, the result for the total energy changes at most by 1.1% 
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10~ 3 
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lO" 4 
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388.23 
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E = Elm, A 


12 


388.08 


14 


388.05 


16 


388.02 


18 


388.00 


20 


387.98 


25 


387.94 


30 


387.90 


40 


387.86 


50 


387.83 



N s =N t 


E = E/m A 


100 


383.35 


120 


384.89 


150 


386.43 


180 


387.46 


200 


387.98 


250 


388.98 


300 


389.52 


500 


390.28 


800 





TABLE I: Energy dependence on numerical parameters. A dash indicates that the relaxation did not 
converge in the predefined number of steps (50000). 



5 10 15 20 25 30 35 40 45 50 

m, R 



FIG. 3: (Color online) Dependence of the total energy (left panel) and the energy per unit length (right 
panel) on the torus radius, at otherwise fixed parameters. 



when varying any of the numerical parameters. Combined with the systematic error from the 
relaxation algorithm (which may not find the absolute minimum), we estimate an overall accuracy 
of about 2% for our classical energy calculations. (This estimate is not indicated by error bars in 
the plots below.) The conclusion is that the numerical parameters in eq. (27) are sufficient to reach 
the continuum limit, at least as long as the particle masses and R are not chosen to have extreme 
values. 

The right panel of figure [3] shows that the straight line (or NO) limit is reached for radii only 
slightly larger than the Compton wave-length of the gauge field. This suggests that energy loss 
from radiating gauge and/or Higgs fluctuations is indeed irrelevant for simulations of networks of 
closed strings that refer to cosmological length scales. 

Next, we want to demonstrate that the torus configuration has a classical instability against 
shrinking into nothing (R — > 0). To see this, we leave the discretization parameters at their values 
in eq. (27) and fix j3 = 1 (and e 2 = ^). 4 In fig. [3J we see that the total energy is a linear function of 
the torus radius R down to very small radii, of the order of a few Compton wave lengths. A better 



4 For very small radii, when opposite sides of the torus start to interact, we need a better resolution in the area 
inside the torus and near the torus core. To do so, we must enlarge the effective r- window and thus [d m in, d ma x] 
while also increasing the lattice extensions to keep the discretization errors small. We have varied the numerical 
parameters in such cases until a stable result was found. 
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FIG. 4: (Color online) Left panel: Higgs profile f(a, r) as a function of a for R = 8.0 and various values 
of r. Right panel: Comparison of the NO Higgs profile from section [TT] with the profile from the torus 
configuration for two values of R. Both charts refer to the model parameters e 2 = \ and j3 = 1. In the right 
panel the variable r refers to r as in eq. |7]) for the NO string and to 2R/(e T — 1) for the closed string. 

signal is obtained by studying the total energy per unit (core) length, E/(2ttR). As can be seen 
from the right panel in fig. [TJ the energy per unit length saturates at a constant asymptotic value 
(which equals the energy per unit length of the NO string listed after eq. (|11|)), down to about 
R < 5. For smaller radii, the interaction between opposite sides of the torus leads to a considerable 
drop in the energy per unit length (cf. density plots below). As a consequence, the total energy E 
of the torus configuration decays more than linearly at small core radii R — > 0, which means that 
it is classically unstable. 

We can also verify that the profiles are predominantly functions of r, with only a very small a- 
modulation. This is shown in the left panel of figure [4] for the Higgs profile f(a, r) (the gauge 
boson behaves similarly). A notable variation with a is only visible at very small r when a ~ 
(i.e. on the vertical axis near the origin). The right panel of fig. [4] demonstrates that the NO 
profiles obtained in the last section agree with the torus profiles in the limit i?> 1, provided that 
we identify the radius r of the NO configuration, eq. ([7]), with the outer distance 2R/(e T — 1), cf. 
eq. p6|). 



Finally, we visualize our configuration with density plots of the Higgs magnitude l^l/m^ ~ f(a, r) 
in a plane perpendicular to the torus core. 5 (The real configuration is found by rotating the 2D 
density plot around the vertical axis.) From the Higgs magnitude, we see the expected profiles for 
large radii, and the overlap in the inner region of the torus for very small radii. The other local 
quantities (magnetic field, energy density etc.) show a qualitatively similar behavior. 

IV. QUANTUM STABILIZATION 



We have seen above that the energy of the torus decreases linearly for a wide range of torus 
core radii R (and even faster than linearly) at very small R), i.e. there is a classical instability of 



5 Here and in the following, the toroidal coordinates are translated back to Cartesian distances for better visualiza- 
tion. 
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the torus configuration against shrinking into nothing. As argued from the uncertainty principle in 
the introduction, we expect that quantum effects will eventually stabilize the torus at a sufficiently 
small radius, so that a stable quantum torus configuration emerges. 



A. Quantum action for the instable mode 



There are various ways to compute the relevant quantum corrections. In situations where the num- 
ber of external degrees of freedom or quantum numbers is large, the leading corrections come from 
the small amplitude fluctuations, which can be efficiently computed using spectral methods [35], at 
least when there is sufficient symmetry to formulate a scattering problem. Formally, this correction 
is of order ff(K). In the present case, however, the instability occurs in a very specific channel, 
associated with the scale of the configuration. This strongly suggests that collective quantum 
fluctuations in this channel may be significant and eventually prevent the configuration from col- 
lapsing. To identify these quantum fluctuations, we promote the torus core radius R to a dynamical 
variable by making it time-dependent, R = R(t). 

More precisely, we consider a time-independent field configuration B/j(x) (not necessarily the 
one that minimizes eq. (19)), where stands for any of the fields in eq. (19). The subscript 



R indicates that the configuration may depend implicitly (through its functional form) on the 
torus radius R. It is important to note that this dependence does not arise from the use of 
toroidal coordinates, which is merely a parameterization of the space point x, but rather describes 
a parameter in the functional form of ®r which is independent of the coordinates used. If we now 
promote R to a dynamical variable by scaling R — > X(t)R = R(t), we obtain a time-dependent field 
configuration 



e(x,t) = e i?w (x), 



(28) 



where x is time- independent. The general form of the action associated with the collective coordi- 
nate R(t) must be 



s = J d 4 xj?[e m ] = J 



dt 



\u(R)R 2 



E(R) 



(29) 



Here, the dot indicates a time derivative and we have explicitly indicated that the functions u and E 
depend on the collective coordinate R(t). In addition they also depend on the model parameters 6 
and the choice of profile functions. For a time-independent radius, the action reproduces the 
energy functional in eq. (19). Moreover, the time derivative R appears quadratically because of 



time-reversal invariance and because the Lagrangian has at most two field derivatives. 

The profile functions depend on the spatial coordinates only via the ratio x/i? that defines the 



toroidal coordinates, cf. eq. (12). We elevate this radius to a time dependent variable. This time 



dependence produces non-zero time derivatives in the action, eq. (29), from which we read off the 



mass function u(R). In addition to the scale dependence of the profile functions the fields may be 
subject to explicit scaling as well. 

The vacuum expectation value of the Higgs fields is fixed at spatial infinity. Thus there is no 
explicit scaling in this case and the parameterization of the time dependent Higgs field becomes, 



(x,t) = fo(x/fl(t)) 



cr t ,T t ) . 



(30) 



6 In units of the gauge mass, u — u/tua and E = E/tua both depend on the gauge coupling only by an overall 
prefactor of 1/e 2 . As explained earlier, this means that the quantum energy scales trivially with e and the optimal 
profiles will depend on the coupling constants only through the mass ratio /3, cf. eq. |9|). 
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10 15 20 25 3D 



FIG. 5: (Color online) Density plots of the Higgs magnitude in a plane perpendicular to the torus. All 
positions are in units of the gauge boson Compton wavelength m^ 1 . In the left panel, the torus has a radius 
R = tuaR = 8, while in the right panel shows results for R — 0.7. Notice the different scales on the axes in 
both cases. 



While <po(cr, t) is a prescribed profile function of toroidal coordinates, at and Tt denote the particular 
toroidal coordinates for the scaled point x/R(t). Since we consider the field </>(x, t) at a given space 
point x, we must treat the arguments of the profile function 4>q as time-dependent variables. As we 
anticipated, this dependence induces a non-zero time derivative that is most conveniently computed 
from eq. (12) for ip = 0, taking into account that dx/dt 



0: 



R 
R 




sincit coshrj 
cosot sinhri 



(31) 



The subscript on the square bracket indicates that the toroidal coordinates are to be taken at the 
time dependent values a = at and r = r^. Then the time derivative of the Higgs fields becomes 



(x,t) 



R 
R 



V<M<r,r) 



where we have defined the dimensionless radial derivative 

d 



V 



d 



sin a cosh r — h cos a sinh t 

da or 



(32) 



(33) 



which does not depend on the radius R. For the gauge field the situation is a bit more complicated 
because there is no boundary condition that disallows an overall scaling by an R dependent function. 
In principle any parameterization would be allowed, because as a solution to the full field equations 
the correct R dependence would be a result. However, we consider only a subset of these equations, 
specified by our ansatz. As for the NO string the field equations require that the covariant derivative 
vanishes as spatial infinity. On the other hand we do not want to abandon our choice g(a, r) — > 1 for 
the boundary condition at |x| — > oo. Since the derivative induces a factor 1/R due to ^ = q^/r) 
we include exactly this factor in the parameterization of the gauge field 



A(x,t) 



1 

W) 



(34) 
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to ensure that the covariant derivative has a uniform scale dependence. This factor additionally 
contributes to the time derivative of the gauge field 



A(x,f) 



R 



(V-l)A (<7,r) 



(35) 



To obtain the Lagrange function, eq. (29), we need to integrate over coordinate space with the 
measure 

d 3 x = R 3 d 3 [x/R] = R 3 d 3 [a t ,T t ,ipt] 



(36) 



where the last expression refers to the integration measure of toroidal coordinates. After spatial 
integration the scale dependence factors as powers R k and the ratio R/R. In summary, the ansatze 
for the field configuration together with their time derivatives are 



$(x,t) 
A(x,t) 



v e 



far), 

n g(oy) 

e Rt](ct,t) e<T ' 



A(x,t) 



v e 



n R 
e R 



■ R 
R 

V- 



inf(a,r) sinu coshr + Vf(a, r) 



g(°y) 



(37) 



Rt](<t,t) 



We have omitted the index on the toroidal coordinates because we have turned them into (dummy) 



integration variables, as in eq. (36). Also, we have not explicitly indicated the implicit dependence 



of the profiles on the radius induced by the stationary conditions. Note that to facilitate the 
numerical calculation, the scaling has been set up such that these profile functions obey exactly 



the boundary conditions discussed in section III B 

The kinetic term A 2 for the gauge boson thus becomes 

R 2 



n\ 2 

e) 2R A r} 2 L 



cos(2<r) + cosh(2r)) g + 4 cosu coshrgVg + 2 CVgY 



(38) 



B. Gaufi' law 



For time-dependent fields, the initial Lagrangian eq. ^ can be split into three pieces: 
& = [\ A 2 + |^| 2 } + {\ (VA ) 2 + e 2 |$| 2 A 2 - A Jo} - {|£>$| 2 + V{$)} 



'stat 



(39) 

The A field only appears in Jz^ aU ss; & n d this contribution gives rise to GauB' law, eq. ([6]), with 
the source 



Jo = ie $) + V A 



(40) 

In the static case, we had Jo = and hence Ao = and Jz? ga uss = 0. This is no longer the case 
when the fields are time-dependent because of the dynamic scaling of the closed string. From <3? in 
eq. (37), the Higgs contribution to the source Jo reads (recall itla = \ 2ev) 



R 



— m' A — f sin a cosh r . 



(41) 



The gauge field contribution to Jo is a bit more involved. Since our torus ansatz in eq. (16) is not 



transverse, we can either take the time derivative of eq. (17) following the rule (37), or take the 
divergence of A in eq. (37). Using the identity [V, V] = —V, the result from both calculations is 



V- A 



n R 
e R? 



V 



(42) 
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From eqs. (41) and (42), we observe Jo ~ R/e and hence eAo ~ R due to the field equation 



(44). The time dependent scaling thus induces a temporal component of the gauge field that is 

the radh 

eAo(x,t) 



proportional to the time derivative of the radius. It is thus suggestive to set 

R 



R 



with a dimensionless profile function q that solves, in principle, 



R 



(-A)+m 2 A f 2 q = J e^ = Q[f,g). 

R 



(43) 



(44) 



Instead of solving eq. (44) directly, it is much more efficient to treat GauB' law on equal footing 



with the remaining field equations. That is, we keep q as an independent field that appears in the 
contribution 



J d 4 xJ? q = J d*x l - {(Vq) 2 + m\ f q 2 ) - qQ = j ' dt± u q (R) R 



(45) 



to the kinetic energy in eq. (29). Since the remaining parts of the action are independent of q, the 
minimization of the action with respect to q will minimize the kinetic term and thus u q , which is 



equivalent to GauB' law eq. (44). The conclusion is that treating /, g, and q as independent fields 



will automatically give the correct contribution to the (quantum) energy from GauB' law, provided 
that q is determined from the minimization of u q (R). 



If we now collect all the pieces from eqs. aS7h, (38) and (45), with Q determined from eqs. (41) 



and (42), the action for the collective coordinate R will indeed be of the form (29). The static 



energy is then given by eq. (19) and the explicit form of the mass term reads 



u(R) 



u(R) 
m A 



2vr 



dr / da sinhr< 



ii 



r ] M[g] + Rg 3 T[f] + 
2ti 



+ Rf 2 v 3 q 2 + l [(d a q) 2 + (d T q) 2 ]-2Rr J 3 JC[f,g)q 



(46) 



The explicit form of the local expressions M, T and JC = Q/m 2 A can be found in appendix [3| To 
discuss the boundary conditions for q(s,t), it is convenient to rewrite eq. (46) in the form 



u(R) 



2tt 

72" 



n 



dt I da sinhW —nM\g] + Rrf (V/) 
2R 



2q^N 
R 2 V 



2) [v' 3 d a (g V )} 



+ 



+ " [(d a q) 2 + {d T q) 2 ] + Rr] 3 n sin a coshr + q ~ f 



R 



(47) 



At very small and very large r, the first line in this equation is suppressed because either the 
derivative of the profiles / and g, or the profiles themselves vanish. The boundary conditions for q 
are thus determined from the balancing of the two terms in the second line of eq. (47). At r — > 0, 
we have / — > 1 by the classical equation of motion, and the last term in eq. (47) dominates. This 
limit thus yields the estimate q 
condition 



-nsincr coshr, which is compatible with the Neuman boundary 
= (48) 



dq(s,t) 
dr 
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also obeyed by the other profiles / and g. At r — > oo, we have / — > by the boundary conditions 



and thus the first term in the second line of eq. (47) dominates, which also yields a Neumann 
boundary condition 



dq(s,t) 



Or 



0. 



(49) 



While this result differs from the Dirichlet condition eq. (22) obeyed be the other profiles, the 
numerical studies below indicate that q(s, oo) when the field equations are obeyed (at the end 
of the (/-relaxation) . 



C. Quantization of the collective coordinate 



Up to now, we have only determined the action (29) due to the unstable mode R(t). To compute 



the quantum correction to the energy due to this mode, we next have to quantize this model. First, 
we note that the mass term u(R) is coordinate-dependent, so we have to introduce a metric factor 
9n = u(R) with g = det(g) = |u(i?)|. The generalized canonical momentum from the classical 



action eq. (29) is 



p = u(R) 



dR 

di 



The coordinate R requires an extra connection (Christoffel symbol, T 



l %u 
29 dR 



upon quantization 



P 



-ih 



— r + -r 

dR 2 



-ih 



d 1 u'(R) 
dR 4 u(R) 



This quantized momentum obeys the canonical commutation relations [R,p] = ih and the Hamilton 
operator of the model (29) becomes 



H 



llA pg 1,2 g ll pg- l/A + E{R) 



p 2 



+ E(R), 



(50) 



where the transformed momentum is p 



pu 



In principle, we would have to find the ground state energy Eq of the Hamiltonian in eq. (50), for 



given background profiles f(a,r) and g(a,r) (with q(a,r) constructed from eq. (H4J)), and then 
minimize Eq rather than eq. ( |19| ). This is very expensive computationally, because we need to 
solve for Eq at every step of the relaxation, which means solving a variational problem for a highly 
non-local functional. As motivated earlier we will follow a simpler approach: From the definition 
of p and p, we find the commutator 8 

[p, R] = -ihu- 1 / 2 , 
and therefore the general uncertainty relation in an arbitrary state tp, 



a^{p) 2 a^(R) 2 > 



h 2 



(51) 



This form of H corresponds to a particular operator ordering which differs from the Weyl ordering. Eq. ( 50 1 is 



instead the standard form in curved space-time, which is usually determined by requiring (i) invariance under 
general coordinate transformations and (ii) reduction to the standard Laplace operator in flat space. 
Here and in the rest of this section, we make h explicit for illustrative purposes. 
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Initially, the Hamiltonian and the wave functions are only defined for R > 0, with a potential 
singularity at R = that is avoided self-consistently. Formally, we can then extend H and ip to 
negative R in such a way that parity is conserved and the ground state ipo is symmetric. As a 
consequence ( ipo | p | i[)q ) = ( | -R | V'o ) = and eq. (51 ) becomes 



< V>o | P 2 | ) > 



fi 2 



4 < | ^ 2 | V'o 



(-00 I u 1/2 1 V>o) 



Next we assume that the ground state ipo is localized such that R is restricted to an effective range 
-R <R<Ro, so that 



(ik\R 2 \ik 



Rn 



and 



(V'o I u~ 1/2 1 Vo) ~ w(-Ro) 



Using this in the previous equation gives 



■00 I P 2 I 00 ! 



fi 2 



u(Ro) Rq 

Then the ground state energy is well approximated by 

h 2 



E = ( Vo I H | 0o > 



2u(R )R 2 



+ U \E{R)\i> 



n 2 



2u{R )R 2 



+ E(Ro) . (52) 



In principle, i?o is the radius of the region where the ground state is localized. For a fixed externally 
prescribed radius R, the ground state will adapt and localize around Rq ~ R. This allows to re- 
interpret eq. (52): For an externally prescribed core radius R, the quantum system is localized 



near R. As with any quantum system, it reacts to the localization with quantum fluctuations that 
are repulsive at very small distances. These fluctuations compete w ith the classical energy E(R), 



which favors R — > according to the instability observed in section 
and the system stabilizes at the minimum of the quantum energy 



IV 



E n 



1 



r, 2 



2 u(R) i?2 



+ E(R) 



These two effects balance 



(53) 



This quantum energy differs from the classical energy by the contribution from the quantization of 
the unstable mode. We assume that this contribution is the dominant quantum effect in stabilizing 
the closed string configuration against shrinking to zero. 



The derivation of equation (53) involved various approximations, but a full quantum mechanical 
treatment is likely to produce something similar — maybe with a numerical factor of order one in 
the quantum correction — because the derived form of the quantum corrected energy is essentially 
dictated by fundamental principles of quantum mechanics. For a first estimate and a proof of 



principle of quantum stabilization, eq. (53) therefore represents a completely adequate starting 
point. 



D. Numerics 



For the numerical minimization of the quantum energy in eq. (53), we employ the relaxation 
technique presented in the last section as follows: We 



1) determine the profiles / and g by minimizing the classical energy E(R) from eq. (19); 
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. E.(R) 



FIG. 6: (Color online) Left panel: Mass factor of the collective coordinate and classical energy, as a function 
of the torus radius. Right panel: Quantum energy vs. classical energy as a function of the torus radius. All 
simulations were performed with the standard parameters e 2 = 1/2 and ft = 1. 



2) compute the auxiliary field q by minimizing the ((/-dependent part of the) collective mass 
function u q in eq. (46) (see appendix [4] for explicit expressions, in analogy to eq. (25) q s t is 
the discretized version form of q(cr, r)), 



5u 
5 q s t 



4"7T 

—iy sinh(r + tA t ) 



X 



st 1st 



St 



0; 



(54) 



3) compute the quantum energy E q (R) from eq. (53) by substituting the profiles from 1) and 2). 



These operations are performed for each fixed radius R separately. Minimization of the total 
quantum energy, E q from eq. (53), would require us to re-adjust the profiles / and g after running 



through steps l)-3). This procedure would then have to be iterated to convergence since modifying 
the profiles / and g causes q to change as well. Alternatively, we could relax all fields /, g and q 
simultaneously using u q as target functional for the minimization of q and E q as target functional 
for / and g. In both cases, the minimization of E q would be performed with the relaxation condition 

h 2 

5E q = ^ 5u + 5E = . (55) 

2u 2 R 2 

Since u is at most quadratic in the profiles / and g, the variation 5u is at most linear. Hence the 



quantum correction in the relaxation condition (55) affects the coefficients H s t, L s t (for /), and 
A s t, B s t (for g), cf. appendix [2] Obviously these corrections are proportional to 1/u 2 and thus 
heavily suppressed since our numerical studies below yield u > 100 for all radii at the minimum, 
and even much larger during relaxation. Therefore the impact of the quantum correction on the 
profiles / and g is a negligible higher order effect and the simple procedure l)-3) listed above is 
essentially equivalent to the complete minimization of E q within our numerical accuracy. 

Figure [6] shows the numerical results for the mass function eq. ( 46 ) for the collective mode using 
the profiles obtained according to steps 1) and 2). For large and moderate radii R, the mass 
is fairly close to the classical energy, and decreases accordingly when R decreases. If this trend 
continued down to very small radii, we would observe quantum stabilization at the expected torus 
size of a few Compton wave lengths, R ~ 1. However, figure [l] shows that the mass function starts 
to saturate 9 at around R ~ 5 at a large value of u ~ 110. As a consequence, the stabilizing radius 



It even slightly increases at very small radii, but this could also be a numerical artifact as the computation becomes 
delicate when R — > 0. 
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R w 0.1 in the right panel of figureJHjis about one order of magnitude smaller than naively expected. 
Although the stabilizing mechanism works in principle, the geometry of the stable configuration is 
characterized by such a small radius that it no longer resembles a torus. It is suggestive that other 
quantum effects may have already become relevant for such configurations. 

A more thorough understanding of the numerical results for u is gained by separating three 



contributions in eq. (46): The first two terms are the gauge and Higgs kinetic energies, respectively. 
The third contribution implements GauB' law. Table [TT] shows some typical results: The gauge 
kinetic energy accounts for about one third of the total mass, while the Higgs kinetic energy and 
GauB' law give very large contributions that almost cancel each other. In total these two pieces 
constitute approximately two thirds of the mass. On the other hand, if we separate contributions 
according to eq. (47) no cancellation of large terms emerges. The reason is that those terms 
are organized as the square of the covariant derivative, l-Do^l 2 . In that combination GauB' law 
guarantees that there are no large canceling contributions to the integral from spatial infinity. If 
we had not solved GauB' law, a non-zero Higgs field contribution to the mass function at spatial 
infinity would not be canceled and the mass function would be ill-defined. This is reflected by the 
(almost linear) increase of the Higgs field component displayed in table |H) 



R = R/rriA 


gauge 


Higgs 


Gaufi 


total 


50.0 


415.0 


227075.0 


-225504.0 


1985.8 


20.0 


167.3 


90108.0 


-89477.2 


798.1 


10.0 


86.0 


44783.1 


-44463.3 


388.7 


8.0 


70.2 


35757.6 


-35499.2 


310.0 


5.0 


47.8 


22260.0 


-22091.4 


216.5 


3.0 


36.2 


13304.9 


-13189.9 


151.2 


2.0 


33.5 


8856.5 


-8759.6 


130.4 


1.0 


32.3 


4435.2 


-4353.3 


114.2 


0.8 


32.3 


3555.9 


-3476.3 


111.9 


0.5 


32.8 


2243.4 


-2164.9 


111.2 


0.3 


31.7 


1323.4 


-1242.5 


112.5 



TABLE II: Contributions to the mass function eq. (46) for the collective coordinate, as a function of the 
torus radius. Data was taken on 350 x 350 lattices at the standard parameters e 2 = 1/2 and (3=1. 



The accurate determination of the mass factor eq. (46) requires a high resolution in the (tiny) 
region between the torus core at R and the origin. In addition, the overall grid extension must 
be large enough to reach asymptotic distances. The numerical parameter d max , which controls 
both the extension towards spatial infinity and the distance from the symmetry axis, is measured 
in units of R. On the other hand, the physical length scale is the Compton wave length, 1/rriA- 
Hence d max has to be very large, especially when R < 1. In practice, values of <i max ~ 200 and grid 
sizes up to 500 x 500 were necessary to obtain stable results. 

For completeness, we have plotted the induced profile function q, that emerges from GauB' law 
for two different radii in figure [7j Deviations from the vacuum value q = extend much further 
out from the core region than the profiles / and g, cf. figure [5j Also, the main structure is located 
above and below the plane of the torus (core line) and q changes sign across this plane. The g-field 



vanishes in the core plane, which is not implemented as a boundary condition, cf. eq. (49), but 
rather a result obtained when the relaxation has converged and thus it is a consequence of GauB' 
law. 

For very small radii, rather large values of d max are required in order to cover the region near 
the vertical axis in figure [7j in which the deviation from the vacuum is sizeable. For illustration 
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2 4 6 8 10 12 



FIG. 7: (Color online) Density plots of the profile function q in a plane perpendicular to the torus. In the 
left panel, the torus has a radius R = 8, while in the right panel shows results for R = 0.7. Notice the 
different scales on the axes in both cases. The profile is positive (negative) in the lower (upper) half-plane. 

purposes d max has been chosen fairly small (d max ~ 10) in figure [7] while in the actual numerical 
calculations it was taken at least an order of magnitude larger. 



V. CONCLUSION 

The Nielsen-Olesen string in scalar electrodynamics is a suitable toy model for cosmic strings 
in non-Abelian gauge theories such as the Standard Model. For this reason we have investigated 
the stability issue of the closed Nielsen-Olesen string. The closed string has been constructed 
by introducing toroidal coordinates. This complicated matters significantly because the profiles 
become functions of two independent variables, and quantitative analyses cannot be pursued in 
from of ordinary differential equations. Therefore we adopted a lattice approach and performed 
variational calculations by relaxation methods. This verified the classical instability of a toroidal 
string as the closed string shrinks to size zero. As in the case of the classical hydrogen atom, 
by radiating energy as fluctuations and simultaneously shrinking, the closed string will eventually 
decay. 

In the context of the cosmological assumption that large scale closed cosmic strings could have 
existed at early times, this analysis suggests that (non- interacting) closed strings would not have 
survived. We have then argued that quantum mechanics in form of the uncertainty principle 
disallows objects that at some stage possessed a finite size squeeze to arbitrarily small extensions. 
Motivated by this observation we have developed a formalism for the quantum stabilization. It 
is based on the uncertainty principle and inspired by its success in reproducing the fundamental 
quantum mechanical properties of the hydrogen atom. 

Technically this approach requires the computation of an inertia function for a time dependent 
radius of the string. This inertia function can be interpreted as a mass for the radius of the torus and 
we have estimated its kinetic energy term in the quantum energy functional. We have then applied a 
variational principle to this functional (appropriately extending the configuration manifold properly 
to implement Gaufi' law) and determined the field configuration from the minimal quantum energy 
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at a prescribed radius. As expected from the uncertainty principle, the kinetic energy part is 
negligible for large radii. More surprisingly, it is also quite small for radii that are of the scale 
of the Compton wave-length of the fundamental fields. Nevertheless we have observed that the 
proposed mechanism of quantum stabilization is indeed operative in the case of the closed Nielsen- 
Olesen string as the quantum energy exhibits a global minimum for a non-zero radius. The proper 
incorporation of Gauss' law proved essential for this mechanism to work. However, the numerical 
results for the radius at at which the torus becomes stable suggest that other mechanisms should 
be relevant as well, since this radius is about an order of magnitude smaller than the length-scale 
inherent in the model. (We do not expect the corrections to the approximations in section IV C to 
compensate for this.) At such small radii the toroidal character of the shrinking configuration has 
disappeared and presumably additional quantum modes (besides the length scale) must be included 
when computing the quantum energy. This conclusion is further supported by the fact that the 
mass associated with the radial mode is quite large. Hence the spectrum of radial excitations will 
be very dense. Presumably these additional quantum fluctuations will be enhanced in the interior 
region of the string. Rather than decaying to infinity, these fluctuations would increase the energy 
of the string thereby stabilizing it at radii larger than those computed here. This suggests that the 
stabilization radius obtained from the application of the uncertainty principle to the radius of the 
closed string is merely a lower limit. Of course, due to the complicated structure of the background 
configuration any quantitative analysis of this scenario will be a formidable task. 
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Appendices 



In these appendices we present details of the numerical approach used in sections III and IV 



1. Lattice discretization 



We start with the discretized description eq. (24), where the profiles turn into real functions 



fst = f(s A s ,r + t A t ) , g st = g(s A s ,r + t A t ) (56) 

defined on the lattice sites. The boundaries of the parameter regions are at t = and t = Nt, and 



at s = and s = N s , respectively, where we impose the following conditions, cf. eqs. (22)-(21): 

• periodic boundary conditions on s: 

fot = fN 3 t , got = 9N a t , t e {0, . . . , N t } (57) 

• Dirichlet conditions at t = Nt: 

f sNt = g sNt =0, se{0,...N s } (58) 

• Neumann conditions at t = 0: 

fso = fsi , 9s0 = 9 si , s e {0, . . . , N s } . (59) 
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Notice that the derivatives in the Neumann conditions are forward derivatives at t = so that no 
extra grid points are required. For integrals over s and t, we use the same prescription combined 
with a left Riemann sum definition, 



y Nt-i 
/ drn(r)^ A t n(r + tA T ) 
{ t=o 

TO 



This ensures that we will never encounter grid points outside the range defined in eq. (24). Next, 
we choose the discretized profiles f s t and g s t such as to minimize the (discretized) target functional 
eq. pb, 



N t -XN a -\ 



t=0 s=0 



In 2 


1 


9s,t+i — 9 st 


\R 


Vst 





+ 



(60) 



+ Rrj st 



fs,t+l — fst 

A, 



+ n 2 f st (l-g st f + 



fs+l,t — fst 

A, 



+U 3 v 3 st(i-f s 2 A, 



where the dimensionless metric factor is given in eq. (20). Once the solutions f s t and g s t are found, 
the original profiles /(x) and g(x) can be restored from eq. (56) and the definition of toroidal 
coordinates, eq. (fl2l). 



2. Relaxation 



For the relaxation condition, we consider an interior point (s, t) at which the profiles are not fixed 
by boundary conditions. At these points we can independently vary the profile functions and 
compute the associated gradients. We start by listing the result for the gauge field: 



SE 2ir 



-1st 



n 2 sinh(ro + tAt 



A st g st - B, 



si 



with 



A 



st 



e ( 1 . it 
R \Vst Vs,t-i 



| +Rf 2 tr]st A 2 



B 



st 



e ( 9s,t+l c 9s,t-l \ A ,2 a 2 
+ it z I + R fst Vst A . 



R V Vst ' Vs,t~i 
For later convenience, we have introduced the abbreviations 
sinh(r + (i - 1) 



it 



sinh(r + tA t ) 



cosh(A t ) — coth(ro + tAt) sinh(A i 



A, 
A t 



(61) 



and A z = A s A t . 

The local minimization of E with respect to g s t yields the updating rule {relaxation condition) 

9st^^, (62) 
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which in view of A s t > is a minimum. For the Higgs profile, f s t, a similar calculation gives 

5E 2vr 



Sf, 



St 



sinh(r + tA t ) R 



^Gst fst + Hst fst — L st 



(63) 



with the coefficients 



G 



St 



P p2„3 A 2 
T-K Vst^ ! 



H st = e (rjst + ^ t -i it) + - (Vst + Vs-i,t) + n 2 i] st (1 - 9st) 2 A 2 - ^ R 2 vlt A 2 , 

L s t = £ (fs,t+l Vst + fs,t-l Vs,t-1 6) + - (fs+l,t Vst + fs-l,t Vs-l,t) ■ 

These values are such that the relaxation condition, fj^ = 0, always has a real solution fjfj!' which 
can be computed by a closed formula into which we substitute the numerically obtained data for 
the coefficients G s t, H st and L st . The updating rule for the Higgs field is then 



f _v f (0) 

Jst r J s t ■ 



(64) 



Since the relaxation conditions (62) and (64) only compute a local minimum with all other field 
variables held fixed, we must proceed iteratively and sweep through the entire lattice many times 
to relax to an overall stable minimum. 



Mass factor for the collective coordinate 



The action for the collective coordinate R, eq. (29), contains a mass factor u that is a local 



functional of the profiles. The three abbreviations in the explicit form eq. ( |46[ ) read 
M = (cos(2a) + cosh(2r)) g 2 + 4 coscr coshr gVg + 2 (Vg) 2 
T = (V/) 2 + (n / sin a cosh r) 2 
K. 



-n 



f 2 sin a coshr+ ^- (V - 2) [rp 3 d a (grj)] . 
R 2 



(65) 



Note that the profile functions / and g depend on the toroidal coordinates a and r. They also 
have an implicit dependence on the torus radius R. 

4. Quantum relaxation condition 



The relaxation condition for the auxiliary field q s t related to Gaufi' law was given in eq. (54) in 
the main text. The corresponding coefficients read explicitly 



X st = Rv 3 st f 2 t A 2 + 



l r 
R 
1 



£ t eVs,t-i + - Vs-i,t + (e + - ] ;/,/ 



Y st = Rvt t lC s tA' + ^ 



1 1 

e Ct Vs,t-i Qs.t-i + - Vs-i,t Qs-i.t + e Vst Qs.t+i + - Vst q s +i,t 
R L e e 



(66) 
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where K, s t is the descretized (lattice) version of K, in eq. (65). 
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